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Abstract We report on new attempt to investigate baryon-baryon interactions in lattice QCD. From the 
Bethe-Salpeter (BS) wave function, we have successfully extracted the nucleon-nucleon (NN) potentials in 
quenched QCD simulations, which reproduce qualitative features of modern AW potentials. The method has 
been extended to obtain the tensor potential as well as the central potential and also applied to the hyperon- 
nucleon (YN) interactions, in both quenched and full QCD. 
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1 Introduction 

In 1935 Yukawa introduced virtual particles, pi- 
ons, to explain the nuclear force^, which bounds 
protons and neutrons inside nuclei. Since then the 
nucleon-nucleon (NN) interaction has been exten- 
sively investigated at low energies both theoretically 
and experimentally. Fig. Q] shows modern NN po- 
tentials, which are characterized by the following 
features' 2 ''^. 
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Fig. 1. Three examples of the modern AW po- 
tential for the 1 So (spin singlet and S'-wave) 
state: Bonn 171 , Reid93 [81 and AV18^. 



At long distances (r > 2 feu ) there exists weak attrac- 
tion, which is well understood and is dominated by 
the one pion exchange, while contributions from the 
exchange of multi-pions and/or heavy mesons such 
as p, oj and a lead to slightly stronger attraction at 
medium distances (1 fm < r < 2 fm). On the other 
hand, at short distances (r < 1 fm), attraction turns 
into repulsion, and it becomes stronger and stronger 
as r gets smaller, forming the strong repulsive core^. 
The repulsive core is essential not only for describ- 
ing the NN scattering data, but also for the stabil- 
ity and saturation of atomic nuclei, for determining 
the maximum mass of neutron stars, and for ignit- 
ing Type II supernova explosions^. Although the 
origin of the repulsive core must be related to the 
quark-gluon structure of nucleons, it remains one of 
the most fundamental problems in nuclear physics for 
a long time^ . It is a great challenge for us to derive 
the nuclear potential including the repulsive core from 
(lattice) QCD. 

In this talk, we first explain our strategy to ex- 
tract AW potentials theoretically from the first prin- 
ciple using lattice QCD and present the recent result 
in quenched QCD simulations. We then apply the 
method to various cases including the energy depen- 
dence of the potential, the quark mass dependence of 
the potential, the tensor potential, the full QCD cal- 
culation and the hyperon-nucleon potential. Since we 
mainly present only the results due to the limitation 
of the length, please see the corresponding references 
for more details. 



Received 31 December 2009 
1) E-mail: saoki@het.ph.tsukuba.ac.jp 



No. X 



Sinya Aoki: Baryon-Baryon Interactions from Lattice QCD 



2 



2 Strategy to extract potentials in 
QCD 

Since a potential is a concept of the non- 
relativistic quantum mechanics, it is non-trivila to de- 
fine it in QCD. To find a reasonable definition of the 
NN potentials in QCD, we first consider the S'-matrix 
below inelastic threshold of the NN scattering. The 
unitarity leads to 

S = e 2li (1) 

where the "phase" S is a hermitian matrix. We next 
introduce the equal time Bethe-Salpeter (BS) wave 
function defined by 



(0\N(x + r,0)N(x,0)\™,E) (2) 



where \2N,E) is a two-nucleon eigenstate in QCD 

-fc 2 , and the N(x) is the 



with energy E = 
gauge-invariant 3-quark operator given by 



N{x) 



c q a {x)q b (x)q c (x). 



(3) 



For large r = |r|, the partial wave I of the BS-wave 
function behaves as 



^ ^ sin(fcr-/7r/2 + ^(fc)) 
kr 



where 5 t (k) is the phase of the S'-matrix for the partial 
wave jl!2>« . (Although S t (k) is a hermitian matrix in 
general, we here consider the case that 6i(k) is just a 
number for simplicity.) The above formula says that 
S t (k) is the scattering phase shift of the scattering 
wave. In other words, the BS wave function defined 
above can be interpreted as the NN scattering wave. 

Based on the above fact, we have proposed the 
following strategy to define and extract the TV TV po- 
tential in QCD f14 ' 15 '. We define a non-local potential 
from the BS wave function as 



[E-H ]tp E ( r ) = J d 3 sU(r,s)<p E ( S ) 



(5) 



where H = — V 2 /(2/Zjv) and fi N = m N /2 is the re- 
duced mass of a two-nucleon system. Since the non- 
local potential is difficult to deal with, we expand it 
in terms of derivatives as U(r,s) = V(r, V)<5 3 (r — s). 
The first few terms are given by^ 

V(r,V) = V c (r) + V T (r)S 12 +V hS L-S 

+ {V D (r),V 2 } + ---, (6) 
3 

S12 = — Oi -0(0-2 -r) -(<Ti- cr 2 ) (7) 

where S 12 is the tensor operator and Oi is the spin 
operator of the i-th nucleon. The central potential 



Vc and the tensor potential V T are the leading local 
terms (without derivatives), and thus can be deter- 
mined from the BS wave function at one energy using 
eq.©. Higer order terms in the above derivative ex- 
pansion can be successively determined from BS wave 
functions at different energy. We however expect at 
low energy that the leading order terms, V c and V T , 
give a good approximation of the potential. We can 
estimate an applicable range of energy for the local 
potential approximation, by calculating physical ob- 
servables such as the scattering phase shift from the 
local potential and comparing them with experimen- 
tal values. 

The first result for the NN potentials based on the 
above strategy has been obtained in quenched lattice 
QCD simulations' l4 J, where the lattice spacing a is 
0.137 fm, the spatial extension L is 4.4 fm and the 
pion mass m n is 529 MeV. In FigfJl the central po- 
tential for the 1 S (spin singlet and L = 0) state and 
the effective central potential for the 3 S*i (spin triplet 
and L — 0) state, obtained at k 2 ~ 0, are plotted as a 
function of r. By comparing Figf2]with FigfTJ we see 
that qualitative features of NN potentials are repro- 
duced. Ref.^ has been selected as one of 21 papers 
in Nature Research Highlights 2007^. 



> 



o 
> 



600 
500 
400 
300 
200 
100 




100 



50 - 



- 



e 
.© 

A & -50 



3§0 



Si — 



$ 4 6 * « * • 



_i_ 



_i_ 



_i_ 



_i_ 



0.0 0.5 1.0 1.5 2.0 
a « * » • » » » * • 



0.0 



0.5 



1.0 
r[fm] 



1.5 



2.0 



Fig. 2. The (effective) central potential for the 
1 So ( 3 Si) state at m-n = 529 MeV in quenched 
QCD. 



3 Recent developments 

In this section we report recent developments of 
lattice QCD calculations for baryon-baryon interac- 
tions, based on the method in the previous section. 
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3.1 Energy dependence 

We first investigate the applicable range of energy 
for the local potential determined at k ~ 0. If terms 
with derivatives such as VlsML-S or {V D (r), V 2 } be- 
comes important, the local potential determined at 
k > is different from the one at k ~ 0. From such 
k dependences of local potentials, in principle, some 
of the terms with derivatives can be determined. In 
Fig El the local potential for the 1 S state obtained 
at k ~ 250 MeV (red, APBC) is compared with the 
one at k ~ MeV (blue, PBC) in quenched QCD at 
a = 0.137 fm and nv = 529 MeV^'^. 
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Fig. 3. Comparison of central potentials in 
the local approximation between the periodic 
boundary condition (PBC, blue) and the anti- 
periodic boundary condition(APBC, red) for 
the 1 5'o state at m w = 529 MeV. 



As can be seen from the figure, the k dependence 
of the local potential turns out to be very small. This 
means that the potential obtained at k ~ in FigJ5] 
well describes physical observables such as the phase 
shift 5 (k) from k ~ to k ~ 250 MeV, in quenched 
QCD at a = 0.137 fm and = 529 MeV. 

3.2 Quark mass dependence 

A quark mass dependence of the NN potential is 
shown in FigfJ] where the central potentials for the 
1 S'o state obtained at k ~ at m, = 380, 529 and 731 
MeV are compared in quenched QCD at a = 0.137 

As quark mass decreases, the repulsion at short 
distance (the repulsive core) get stronger while the 
attraction at intermediate distance ( 0.6 fm ~ 1.2 
fm) becomes also stronger. 
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Fig. 4. The central potential for the 1 So state 
at =380 MeV (red), 529 MeV (green) and 
731 MeV (blue) in quenched QCD at a = 0.137 
fm. 



3.3 Tensor potential 

The tensor operator S 12 mixes the 3 S 1 ( J = S = 1 
and L = 0) state with the 3 D 1 ( J = S = 1 and L = 0) 
state. Using this property, we can determine the ten- 
sor potential as follows. For the J = S = 1 state, the 
local potential approximation leads to 

(E-H a )<p B ( r ) = [V c (r)+V T {r)S 1% ]y E {r), (8) 

which, by the projection P to the L = state and the 
projection Q to the L = 2 state, is decomposed into 



Pip E PS 12 <p E 
Qifi E QS 12 tp E 





= (E-H 



The above equation can be easily solved as 

PlPe PS 12 ip E \ 
Qip E QS 12 (p E J 



(9) 




x (E-H ) 



(10) 



In Figl5l the central potential V c (r) and the ten- 
sor potential V T (r) for the spin-triplet state are 
plotted^, together with the effective central poten- 
tial V,S ff (r) for the 3 S'i in FigJ^l which corresponds to 
(E—H )Pip E J (P<Pe) in the above notation. These po- 
tentials are calculated in quenched QCD at a = 0.137 
fm and = 529 MeV. 
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3.4 Full QCD calculations 
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Results for the NN potentials so far have been 
calculated in quenched QCD. In Ref.^, prelimi- 
nary results in full QCD calculations have been re- 
ported, based on gauge configurations generated by 
the PACS-CS Collaboration in 2+1 flavor QCD at 
a = 0.09 fm and L = 2.9 fm™. Hadron spectra ob- 
tained from these configurations are shown in FigjTJ 
Agreements between lattice QCD predictions and ex- 
perimental values are quiet good. The (effective) cen- 
tral NN potentials on these configurations are given 
in Fig|5]for the 1 S state (red) and the 3 S 1 state (blue) 
at TO, = 702 MeV. 



Fig. 5. The central potential(blue) and the ten- 
sor potential (red), together with the effec- 
tive central potential (red), for the spin-triplet 
state at = 529 MeV. They are obtained 
from k ~ states in quenched QCD at a — 
0.137 fm. 

We first notice that the tensor potential V T has no 
strong repulsive core, in contrast to the central po- 
tential. This feature is consistent with the previous 
phenomenological estimate' 21 -^. Although the tensor 
potential is comparable in the magnitude with the 
central potential, the difference between the central 
and effective central potentials, which is caused by the 
second order perturbation of Vt, is very small at this 
quark mass. A quark mass dependence of the tensor 
potential is given in Fig(51 where the tensor poten- 
tial is plotted at = 380,529 and 731 MeV^. The 
tensor potential becomes stronger as the quark mass 
decreases. 
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Fig. 7. Light hadron spectra extrapolated to 
the physical point using m n , vhk and mn as 
input. Horizontal bars denote the experimen- 
tal values. 
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Fig. 6. Quark mass dependence of the tensor 
potential. 
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Fig. 8. The (effective) central potential for the 
1 So state (red) and the 3 Si state (blue) at 
= 702 MeV in 2+1 flavor QCD at a = 0.09 

fm. 
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We observe that the repulsive core in both states 
is much larger in magnitude than the corresponding 
one in quenched QCD in Figj2j though the lattice 
spacing ( 0.09 fm vs. 0.137 fm) and the pion mass ( 
529 MeV vs. 702 MeV ) are different. A reason for 
the difference of the repulsive core between full and 
quenched QCD is now under investigation. 

3.5 Hyperon-Nucleon interactions 

A hyperon is a baryon which contains at least 
one strange quark. Contrary to the case of NN 
interactions, hyperon-nucleon (YN) or hyperon- 
hyperon(W) interactions can not be precisely de- 
termined, since the scattering experiments are either 
difficult or impossible due to the short life times of hy- 
perons. Our approach therefore may open a new pos- 
sibility to determine them theoretically from QCD. 
In this direction, potentials between a S° (hyperon 
with strangeness —2 ) and a proton has already been 
calculated in quenched QCD [24 '. 
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Fig. 9. The central potential for the AN ^So) 
state in 2+1 full QCD as a function of r at 
m w ~400 MeV (red) and 700 MeV (green). 







50 - 






40 - 






30 - 
20 - 




> 

5 


10 - 

- 






-10 - 
-20 - 






-30 - 
-40 - 
-50 - 




r(fm) 

Fig. 10. The central and tensor potentials for 
the AN ( 3 Si - 3 Di) state in 2+1 full QCD at 
m w ~ 400 MeV (red and blue) and 700 MeV 
(green and magenta). 

Recently calculations have been extended to the 
potentials between the A( hyperon with strangeness 



-1 ) and N in both quenched and full QCD. FigU 
and FigHU] show the AV potentials as a function 
of r obtained from 2+1 flavor QCD calculation' 25 ' 
based on the PACS-CS gauge configurations. The 
central potential for the 1 S state is given in Figj9l 
while the central and the tensor potentials for the 
A S 1 — 3 D 1 state are given in FigJTU], with highlighting 
the short distance (medium to long distance) region 
in the left (right) panel. These figures contain results 
at m n ~ 400 and 700 MeV. 

As can be seen from both figures, the attractive 
well of the central potential moves to outer region as 
the quark mass decreases while the depth of these at- 
tractive pockets do not change so much. The present 
results show that the tensor force is weaker while the 
spin dependence is stronger than the JVJV case^ . As 
in the NN case, the hight of the repulsive core is much 
larger than the quenched case^ and it increases as 
the quark mass decreases. 

4 Conclusion 

We present recent results on baron-baryon inter- 
actions obtained from lattice QCD simulations. In 
our strategy, baryon-baryon(W, YN and YY) po- 
tentials are extracted from the BS wave functions. 
The first result for the NN (effective) central poten- 
tials in quenched QCD shows good "shape": Quali- 
tative features of the NN potential have been repro- 
duced in quenched QCD, and the energy dependence 
of the potentials is weak at low energy. The method 
has been successfully extended to the tensor poten- 
tial and the AV potentials in both quenched and full 
QCD. 

One of the ultimate goal in our approach is to 
calculate baryon-baryon potentials in full QCD at 

= 140 MeV. In such calculations one can inves- 
tigate, for example, a relation between the deuteron 
binding and the tensor force. As other directions of 
our approach, it is important to extract the 3-body 
force^ from QCD and to understand the origin of 
the repulsive core theoretically^ . 

/ would like to thank members of HAL QCD Col- 
laboration for providing me the latest results and 
useful discussions. This work is supported in part 
by Grants-in-Aid of the Ministry of Education, Cul- 
ture, Sports, Science and Technology of Japan (Nos. 
20340047, 20105001, 20105003). 
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